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Abstract. We investigate the capability of the UBVRIJHK Yihotometric sys- 
tem to quantify star clusters in terms of age, metallicity and color excess by 

their integrated photometry in the framework of PEGASE single stellar popu- 
lation (SSP) models. The age-metallicity-extinction degeneracy was analyzed 
for various parameter combinations, assuming different levels of photometric 
accuracy. We conclude, that most of the parameter degeneracies, typical to 
the UBVRI photometric system, are broken in the case when the photome- 
try data are supplemented with at least one infrared magnitude of the JHK 
passbands, with an accuracy better than '^0.05 mag. The presented analysis 
with no preassumptions on the distribution of photometric errors of star cluster 
models, provides estimate of the intrinsic capability of any photometric system 
to determine star cluster parameters from integrated photometry. 

Key words: techniques: photometric - methods: data analysis - galaxies: 
star clusters 

1. INTRODUCTION 

Broad-band photometric data of extragalactic star clusters are applicable to 
derive their evolutionary parameters via comparison with single stellar population 
(SSP) models. However, this procedure is restricted by the fact that in some pa- 
rameter domains strong age-metallicity-extinction degeneracies (see, e.g., Worthey 
1994) remain unsolved. For a study of broad-band color indices, most sensitive 
to various star and cluster parameters, see, e.g., Jordi et al. (2006) and Li et al. 
(2007). 

It is known (Anders et al. 2004; Kaviraj et al. 2007) that infrared and ultravi- 
olet passbands are helpful for breaking the age-metallicity-extinction degeneracies. 
In this respect infrared observations have the priority: in some atmospheric win- 
dows they are accessible for ground-based telescopes, while ultraviolet photometry 
shorter than 300 nm is possible only from space. In some extragalactic studies of 
star clusters published to date various combinations of optical and/or near-infrared 
passbands have been used (e.g., Kodaira et al. 2004; Hempel & Kissler-Patig 2004; 
Fan et al. 2006; Hempel et al. 2007; Narbutis et al. 2008; Pessev et al. 2008). 

However, in comparison to photometry in the optical range, ground-based in- 
frared photometry usually has larger errors due to variable humidity, and this 
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requires a careful calibration (e.g., Kodaira et al. 1999; Kidger et al. 2006). The 
reliability of secondary standards used for photometric calibration of wide field 
images, should be verified by several independent sources, if they are available 
(e.g., Narbutis, Stonkute & Vansevicius 2006). 

In the present paper we continue to investigate a possibility of determining 
star cluster parameters (age, metallicity and color excess) by comparison of their 
integrated color indices with the SSP models computed with the PEGASE (v. 2.0; 
Fioc & Rocca-Volmerange 1997) code package. In our previous study (Narbutis 
et al. 2007b, hereafter Paper I) it was found that the UBVRI system enables us 
to estimate cluster parameters over a wide range of their values, when the overall 
accuracy of color indices is better than ~0.03 mag. In the following we discuss 
how the adding of the JHK passbands affects the accuracy of cluster parameter 
(age, metallicity and color excess) determination. We analyze degeneracies at 
various accuracy levels of photometry for the same values of cluster parameters as 
in Paper I. 

In the similar study Anders et al. (2004) have used the so-called AnalySED 
method for determining cluster parameters with a different approach to the degen- 
eracy problem. In Section 2 we discuss the main differences between the Anders et 
al. and our methods and note, that the parameter degeneracy analysis presented 
in this study is based on minimum assumptions. 

2. THE METHOD 

The analysis method used in this study is similar to that of Paper I. The 
SSP models were computed with the PEGASE program package, applying its 
default options and the universal initial mass function (UIMF; Kroupa 2001). 
The integrated color indices in all pass-bands in respect to the t/-band of SSP 
models were reddened by taking into account the dependence of color-excess ratio 
(e.g., Ejj^B / Eb-v) on color index B-V of SSP model and assumed color excess, 
Eb-v, applying the standard extinction law (Cardelli et al. 1989). The three- 
parameter (3-D) SSP model grid of 5 x lO'^ models was constructed at the 
following nodes: (i) 76 age, t, values from 1 Myr to 20 Gyr with a constant step of 
log(t/Myr) = 0.050 (ii) 31 metalhcitjEI, [M/H], values from -2.3 to +0.7, with a 
step of 0.1 dex; (iii) 201 color excess, Eb-v, values from 0.0 to 2.0, with a step of 



The procedure of determination of cluster parameters {t, [M/H], Eb-v) was 
implemented as a C-l — I- code in the data analysis and the graphing software package 
'Origin' (OriginLab Corporation). It is based on a similar technique developed for 
star quantification by Vansevicius & Bridzius (1994), i.e., the comparison of the 
observed color indices of a star cluster with color indices of the SSP from the model 
grid. For this purpose we use the quantification quality criterion, 6, calculated by 
the formula: 



where CI° stands for the "observed" color indices U-V, B-V, V-R, V-I, V-J, 



0.01. 




(1) 



^ For the ages < 12 Myr the step was equal to 1 Myr. The SSP models for the ages larger 
than the oldest known globular clusters are used to avoid marginal effects. 
^ [M/H] was computed applying the approximation [M/H] = log(Z/ZQ). 
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V-H and V-K; C/™°'^ - the corresponding color indices of the SSP models from 
the grid; Wi - the weights for the observed color indices. 

We investigate the possibility to determine cluster parameters using the UB- 
VRIJHK photometric system for 54 models taken from the SSP model grid as 
"observed" objects with the following parameters: t = 0.02, 0.05, 0.1, 0.2, 0.5, 1, 

2, 5, 10 Gyr, [M/H] = 0.0, -0.7, -1.7 and Eb-v = 0.1, 1.0. Their color indices 
are denoted by Clf^^, see Eq. 1. We have assigned the weights of = 1 for U — V, 
B -V,V - R andV - I color indices, and = 1/9 for V-J, V-H and V-K of 
the "observed" star clusters. Such weighting scheme assumes that the accuracy of 
the infrared Jffif observations is ~ 3 times lower than that of the optical UBVRI 
- a typical condition in photometry of extragalactic star clusters. 

Anders et al. (2004) have extensively analyzed the possibilities of their Anal- 
ySED method, which has a different approach to the degeneracy problem, com- 
paring to the method presented in this paper. First, AnalySED is based on the 
comparison of observed magnitudes with the model magnitudes, while for the same 
aim we are using color indices. This allows us to avoid fitting of an additional free 
parameter - the cluster mass. Second, Anders et al. (2004) define the "observed" 
magnitudes of model clusters by assigning to them random errors. However, we 
did not apply "observation errors" to color indices, since it is reasonable to as- 
simie that the quantification quality criterion, 5, reflects the effect of photometric 
accuracy sufficiently well. In general, 6 represents the lowest possible average pho- 
tometric error of all color indices, which in the case of real clusters can be larger, 
mostly due to the influence of crowding by contaminating fleld stars (e.g., Nar- 
butis et al. 2007a). Therefore, the parameter degeneracy analysis presented in 
this study is based on minimum assumptions. 

We have used several S^ax values as the upper threshold levels of photometric 
errors. The cluster parameters were determined independently at each accuracy 
level by averaging the corresponding parameter values of SSP models at the grid 
nodes, which have 5 < (5max- The weights for the calculation of the parameter 
averages and the standard deviations were assigned 1 and 10~^/(5^ for the nodes 
with S < 0.01 and > 0.01 mag, respectively. 

Furthermore, for the parameter averaging we used only the nodes, which reside 
in a single continuous 'island' around the true cluster position in the 3-D parameter 
space of the SSP model grid. Boundaries of 'islands' at each 6max level were 
determined automatically by the clustering procedure, which finds discontinuities 
in the parameter space, starting from the true position of the cluster. Such a 
procedure excludes nodes, which are located in the secondary 6 minima, arising 
due to the age-metallicity-extinction degeneracies in the 3-D parameter space. 

3. RESULTS AND DISCUSSION 

The results of parameter determinations of clusters are provided in Figures 1 
and 2 for color excess values Eb-v = 0.1 and 1.0, respectively. In each panel 
the differences of parameters (determined minus true) are shown in groups of five 
filled circles, the middle circle is positioned at the true age of the cluster, i.e., 
the circles indicate parameters determined at (5,nax threshold values from 0.01 to 
0.05 mag, with a step of 0.01 mag, plotted from left to right. Open and black circles 
correspond to the UBVRI and f/SVS/Jifii'passband set cases, respectively. Error- 
bars indicate standard deviations of the determined parameters and characterize 
the integrated 'size' of the 6 < Smax 'island' in the corresponding parameter space. 




0.1 1 10 0.1 1 10 0.1 1 10 

t [Gyr] f IGyr] t [Gyr] 

Fig. 1. Star cluster parameter determination results for star clusters with the following parameters: age, t = 0.02, 0.05, 0.1, 0.2, 
0.5, 1, 2, 5 10 Gyr; metallicity, [M/H] = 0.0, —0.7, —1.7 (left, middle and right panels); color excess, Eb-v = 0.1. Differences of the 
parameters (determined minus true) are shown in groups of five filled circles, the middle of which is positioned at the true age of a star 
cluster. Within a group, from left to right, the circles correspond to the quantification results derived for (5max from 0.01 to 0.05 mag, 
with a step of 0.01 mag. Open and black circles correspond to UBVRI and UBVRIJHK passhand set cases, respectively. Error-bars 
indicate the standard deviations of the derived parameters. 
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Fig. 3. The distribution of tlie quantification quality criterion, 5, in projection onto 
the Eb-v and [M/H] vs. log (t/Myr) pianos for the clusters with the parameters: Eb-v ~ 
0.1, [M/H] = 0.0, t = 0.2 and 2 Gyr (left and right panels); Eb-v = 0.1, [M/H] = -1.7, 
f = 0.5 Gyr (middle panel) for the case of UBVRI passband set. The true positions 
of cluster parameters are marked by a white "plus" symbol. Black, gray and light-gray 
shaded areas correspond to (5max = 0.01, 0.03 and 0.05 mag, respectively. 



It is clearly visible in both Figures 1 and 2, that additional usage of JHK 
passbands significantly improves the precision of cluster parameter determination 
with respect to the i7i?V72/ passbands alone. Especially good results are in the case 
of the small color excess, Eb-v = 0.1, value. Note, that strong age-metallicity 
degeneracies at ages t = 0.1, 0.2 and 10 Gyr, which are present in the UBVRI case, 
completely disappear in the UBVRIJHK case. A similar, but a slightly weaker 
effect is also visible for models of ages t = 2 and 5 Gyr. However, if objects are 
highly reddened, Eb-v ~ 10, there are still notable parameter degeneracies for 
some cases. The strongest age-extinction degeneracy is for clusters with: t = 1 Gyr 
and [M/H] = 0.0; t = 2 Gyr and [M/H] = -0.7; t = 5 Gyr and [M/H] -1.7. 
Actually, the quantification results of the last case are much more accurate, when 
SSP models older than 15 Gyr are not included in the calculation of the average. 

The impact of adding the Jf/X passbands to the C/SFi?/ system on the accuracy 
of the quantification of cluster parameters is illustrated in Figures 3-10 for the 
clusters of ages t = 0.2, 1, 2, 5, 10 Gyr for solar metallicity, and t = 0.5 Gyr for 
[M/H] = —1.7. The figures display cluster parameter quantification maps at the 
threshold levels of (^max = 0.01, 0.03 and 0.05 mag. 

In Figures 3 and 5 for the UBVRI set, the left and the central panels look very 
similar, although the left panel displays quantification map of a cluster with t = 
200 Myr and [M/H] = 0.0, while the central panel - for a cluster with t = 500 Myr 
and [M/H] = —1.7. This means that due to age-metallicity degeneracy for both 
clusters we obtain similar and wrong parameters. However, the situation changes, 
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Fig. 4. The same as in Figure 3, but for the case of UBVRIJHK passband set. 
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Fig. 5. The same as in Figure 3 for the case of ?7B Fil/ passband set, but for Eb-v 
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Fig. 6. The same as in Figure 5, but for the case of UB VRIJHK paashand set. 
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Fig. 7. The same as in Figure 3 for the case of UBVRI passband set, but for the star 
clusters with the following parameters: Eb-v = 0.1, [M/H] = 0.0, t = 1, 5 and 10 Gyr. 
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Fig. 8. The same as in Figure 7, but for the case of ?7S V7?/Jffiir passband set. 
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Fig. 9. The same as in Figure 7 for the case of ?7B Fil/ passband set, but for Eb-v = 
1.0. 
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Fig. 10. The same as in Figure 9, but for the case of UB VRIJHK passhand set. 



when the JflX magnitudes are added. Figures 4 and 6 display quantification maps 
of the same clusters, but for the UBVRIJHK passhand set. Note, that degeneracies 
disappear, even at Jmax = 0.05 mag level. 

An interesting example is the cluster with the parameters t = 2 Gyr and 
[M/H] = 0.0. Figures 1 and 2 show that the addition of the JHK magnitudes 
improves the age determination, but Eb-v determination becomes of lower accu- 
racy, comparing to the UBVRI set alone. Rightmost panels in Figures 3 6 clarify 
this situation. In the UBVRI case, strong agc-nictallicity degeneracy takes place: 
note the additional 'tail' of suitable SSP models in the Eb-v vs. t plane towards 
older ages and lower Eb~v values. This helps to compensate the influence of SSP 
models at higher Eb-v values, and the quantification produces a 'correct' redden- 
ing. In the case of UBVRIJHK, no such 'tail' is seen, therefore the determined 
Eb-v value is shifted from its true position. 

Figures 7-10 display three similar cases of parameter degeneracies for clusters 
of solar metallicity and ages t = 1, 5 and 10 Gyr. Only the 10 Gyr age sample 
is sensitive enough to the addition of the JHK passbands, when parameter degen- 
eracies become broken even at Jmax = 0.05 mag level. On the contrary, for the 
ages of 1 Gyr (except for the case of Eb-v = 0.1; see Figures 7 and 8) and 5 Gyr, 
parameter degeneracies exist even for (5niax = 0.01 and 0.03 mag. respectively. 
This implies, that the overall accuracy of cluster color indices must be better than 
0.03 mag, which is rather difficult to achieve in infrared photometry. 

We also investigated separately the influence of each infrared passband to the 
parameter determination precision with the purpose to find a minimum passband 
set, which would be sufficient for the quantification of clusters. The result is 
that all three color indices {V-J, V-H, V-K), added to the t/SFii/ photometric 
system separately, reduce the quantification errors, but V-K is most important. 
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Fig. 11. Comparison of star cluster parameter quantification results for the case of color excess Eb-v = 1-0 (see Figure 1 for star 
cluster details). Black filled circles correspond to the case of (7i? VTJ/Ji/A' passband set, open circles - to the UBVRIK. The error-bars 
show the sensitivity of determined cluster parameters to the artificial shift of the V-K color by ±0.03 mag. 
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Figure 11 displays the quantification results for clusters using the UBVRIJHK 
(filled circles) and the UBVRIK (open circles) passband sets, both pictures are 
quite similar. Thus the K passband helps to solve majority of the parameter 
degeneracy problems. The error-bars in the UBVRIK case show the sensitivity of 
the parameters to the change of V-K by ±0.03 mag. 

Recently an updated version of stellar isochrones, compared to those used in 
PEGASE, has been released by the Padova group, which includes a more accu- 
rate treatment of thermally pulsating AGB stars (Marigo et al. 2008). It has 
been shown that the interpretation of integrated colors of unresolved galaxies, 
by applying galaxy populations synthesis method and the new isochrone set, can 
be significantly altered due to improved models of AGB stars (e.g., Tonini et al. 
2008). In the case of single stellar population analysis of star clusters one can 
expect similar changes. Our method of the parameter degeneracy analysis, which 
is much simpler than the methods described in the literature (see, e.g., de Grijs 
et al. 2005), can be used to study the effects AGB stars on cluster colors by 
intercomparing parameter degeneracies in various SSP model frameworks. 

4. CONCLUSIONS 

We conclude that additional photometric information from the JHK passbands 
can significantly improve the accuracy of the determination of cluster parameters 
based on the PEGASE SSP models, in comparison with the results when photom- 
etry only the ?75Fi?/ passbands is available. Even one additional if passband can 
improve significantly the capability of the f/BVT?/ photometric system to eliminate 
age-metallicity and age-extinction degeneracies in the majority of the investigated 
cluster models, when the overall accuracy of color indices is better than ~0.05 mag. 

Note, however, that this condition of photometric accuracy is broken for young 
low-mass star clusters, where the stochastic effects, arising due to a few bright 
stars, dominate (e.g., Cervino & Luridiana 2004; Deveikis et al. 2008). This 
implies, that even in case of ideal photometric calibrations, the uncertainty of 
cluster colors due to stochastic effects limits the applicability of SSP model fitting 
to derive evolutionary parameters of low-mass star clusters. 
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